1
Limitaciones de los Polinomios Característicos
MATH007Lesson 9
00:00
Aunque el polinomio característico $p(\lambda) = \det(A - \lambda I)$ es la base teórica para definir los valores propios, es numéricamente "mal condicionado" y computacionalmente ineficiente para sistemas de alta dimensión. En aplicaciones prácticas —como resolver el sistema de Sturm-Liouville para la propagación de ondas— la sensibilidad de las raíces del polinomio a perturbaciones en sus coeficientes hace que la expansión directa sea una opción secundaria.

De Ondas Continuas a Matrices Discretas

La vibración de una cuerda o membrana está gobernada por la ecuación de onda:

$$\rho(x) \frac{\partial^2 v}{\partial t^2}(x, t) = \frac{\partial}{\partial x} \left[ p(x) \frac{\partial v}{\partial x}(x, t) \right]$$

Para hallar la solución $v(x, t) = \sum_{k=0}^{\infty} c_k u_k(x) \cos \sqrt{\lambda_k}(t - t_0)$, debemos resolver el sistema de Sturm-Liouville:

$$\frac{d}{dx} \left[ p(x) \frac{du_k}{dx}(x) \right] + \lambda_k \rho(x) u_k(x) = 0$$

Complejidad de la Discretización

La discretización del operador da lugar a ecuaciones matriciales como $Aw = -0.04 \frac{\rho}{p} \lambda w$. Para una matriz tridiagonal $4 \times 4$, $p(\lambda)$ es manejable. Sin embargo, a medida que se refinan la malla ($n$ aumenta), nos encontramos con dos obstáculos:

  • Límite de Abel-Ruffini: No existe una solución algebraica para las raíces de los polinomios donde $n \ge 5$.
  • Sensibilidad al Redondeo: En sistemas de alta dimensión, un cambio en la cifra decimal $10^{-10}$ de un solo elemento puede desplazar los valores propios por órdenes de magnitud (fenómeno del polinomio de Wilkinson).

Necesidad Numérica y Bibliotecas Profesionales

Las bibliotecas numéricas profesionales (IMSL, NAG) evitan los polinomios característicos directos. En su lugar, utilizan procedimientos iterativos para aproximaciones:

  • Biblioteca IMSL: Utiliza mínimos cuadrados lineales, splines cúbicos y transformadas rápidas de Fourier.
  • Biblioteca NAG: Emplea aproximación polinómica por mínimos cuadrados y ajustes en sentido $l_1/l_{\infty}$.

Al aproximar valores propios para el sistema $\lambda_i = 1 + 4\alpha\left(\sin \frac{\pi i}{2m}\right)^2$, confiamos en los mínimos cuadrados discretos y en el descubrimiento iterativo en lugar de encontrar raíces.

🎯 Herramienta Teórica frente a Riesgo Numérico
El polinomio característico $p(\lambda) = \det(A - \lambda I)$ es vital para las demostraciones, pero peligroso para los cálculos. Los problemas prácticos de valores propios en física se resuelven mediante transformaciones iterativas (como QR) que preservan la estabilidad.